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An analysis is made of current heat flow data and thermal models of lunar evolution 
which satisfy the diverse information that has accumulated on internal processes. 


Direct measurements of the heat flux were 
made in drill holes at the Hadley Rille 
(Apollo 15) and Taurus-Littrow (Apollo 17) 
regions located on the margins of Mare Sere- 
nitatis and Mare Imbrium, which apparently 
have mascons. These measurements resulted 
in comparatively high values of the heat flux, 

3.1 X 10-8 w/cm2 ( 0.74 X 10-« 5^1 ^ 

\ cm^sec / 

(Apollo 15) and 

2.8 X 10-8 W/cm* ( 0.67 X 10-« — \ 

\ cm-sec / 

(Apollo 17), which differ from the wide 
spread earlier concept of a much lower heat 
flow from a Moon of chondrite composition. 
Only a small number of earlier works pre- 
dicted a high heat flow. Those were the data 
on the thermal emission of the Moon in the 
microwave region (ref. 1) and the results 
of certain calculations of the thermal history 
of the Moon (ref. 2). Current interpretation 
of the high heat flow requires a high average 
concentration of uranium for the Moon — 
60 ppb — as well as an initial surface tem- 
perature sufficient to melt the upper several 
hundred kilometers, including the lunar crust 
and lithosphere. 

The crust of the Moon is a discrete layer 
several-dozen-kilometers thick, bears traces 
of magmatic fractionation that occurred 
during the period of intensive melting of the 
uppermost layers of the Moon (hundreds of 
kilometers deep), and has been reworked by 
the impacts of meteorites. 


Among the sources of initial melting are 
the following : accretional heating during ac- 
cumulation of the Moon, possible electrical 
heating by an early strong solar wind, tidal 
friction when the Moon was close to the 
Earth, and the effect of short-lived isotopes. 
All the sources are problematical. The ther- 
mal effect due to collisions of protomoons 
during accretion appears to be adequate to 
melt the rocks in the area of impact, but not 
at great depth (refs. 3 and 4). 

At this stage of investigation of the 
thermal state of the Moon, the following 
problems are the most significant. 

What was the cause of the rapid early 
melting and differentiation? 

To what extent can heat fluxes obtained 
close to the margin of a mascon be repre- 
sented as global values? To what extent can 
an average, presently uniform concentration 
of uranium be responsible for the surface 
heat flux that undoubtedly results from ir- 
regular distributions of uranium and thorium 
with depth, due to the rapid fractionation 
during the accretionary stage? To what ex- 
tent can coarse numerical modeling of the 
thermal history of the Moon, where the first 
step in the calculation grid is 20 km in size 
(refs. 5 and 6) or where even the smallest 
is 5 km (ref. 7), reflect the fine structure of 
the uniquely low thermal conductivity of the 
upper few meters? 

Before the Luna and Apollo flights to the 
Moon, a lunar model of chondrite composition 
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was widely used (refs. 2 and 8 through 13). 
A drastic reexamination of these models was 
required after study of the chronology of the 
lunar samples and determination that the 
lunar K/U ratio was 2000 instead of 80 000 
as for chondrites and 10 000 as for the Earth. 
The low lunar K/U ratios are in agreement 
with the overall depletion of the Moon in 
volatile elements. 

A comparative study of the energy balance 
of the Earth and the Moon by Lyubimova 
(ref. 2) showed that the integral heat loss 
Qa for the lifetime of the Moon, the integral 
heat generation Ha, and the total heat content 
at the melting point of the planet .\<j are in the 
relationship Aa< Hq ^ Qa for the Moon. This 
means that the energy losses from the surface 
of the Moon predominate (in distinction from 
the Earth) over internal heat generation and 
heat content for the state of complete melt- 
ing, i.e., the Moon is cooling down. The heat 
flux at the surface is assumed to be close 
to the upper limit established by lunar sur- 
face radio emission methods (refs. 14 and 
15). According to thermal history calcula- 
tions, the integral heat flux of the Moon 
should amount to Q = 1.9 X 10®* ergs, while 
the heat content of a completely molten Moon 
is A<j = 1.6 X 10®* ergs (ref. 2). It was con- 
cluded that even in the case of a quite low 
concentration of radioactive elements, in 
agreement with the chondritic model, the 
energy is completely sufficient to convert the 
lunar interior to the molten state. 

Direct measurements of the heat flux dur- 
ing Apollo 15 and Apollo 17 permit one to 
return to discussion of the energetic state of 
the Moon. They also make it possible, to- 
gether with radioastronomers, to refine mod- 
els of the structure of the uppermost layers 
of the Moon (refs. 16 and 17). Serious limi- 
tations on the thermal history of the Moon 
were introduced as a result of a seismic ex- 
periment on the Moon (ref. 18). They in- 
dicate the necessity for a high temperature 
of the lunar interior and the existence of an 
asthenosphere below 1000 km. On the other 
hand, the internal temperature derived from 
the electrical conductivity profile indicates 
that the temperature of the lunar interior 


should be less than the solidus temperature 
(ref. 19). The existence of mascons also re- 
quires this because the thickness of the litho- 
sphere must be sufficient to sustain a load 
corresponding to the mascons. Correspond- 
ingly, the viscosity of the lunar interior is 
assumed to be very high, up to 1(F® poise 
(ref. 20). The chronology of lunar volcanic 
activity shows that all episodes of volcanic 
activity and excavation of the mare must 
have ended 3 b.y. ago (refs. 21 and 22). The 
magnetic history of the Moon indicates that 
there was an old magnetic field which is re- 
sponsible for the remanent magnetism of the 
Moon (ref. 23). This again entails the con- 
clusion that there was intense melting in the 
initial stages of evolution of the Moon. In 
this manner, the factual materials returned 
by the Luna and Apollo missions indicate 
that the Moon evolved quickly. 

In this article, we examine the parameters 
of lunar thermal models, calculation tech- 
niques that consider melting and effective 
convection, and the expected regional correc- 
tions to the heat flux. 


Evolution of a Homogeneous 
Moon 


In a number of the most recent works, the 
heat flux from the lunar interior is compared 
with an average uniform concentration of 
uranium (refs. 5, 6, and 23) in a homogen- 
eous Moon. 

A strict analytical solution for a body of 
spherical shape (the Moon in particular) 
with uniform distribution of heat sources, 
thermal diffusivity k, density p, and heat ca- 
pacity C, was given by Lyubimova (refs. 2, 
24, and 25) , in the form: 


T, 


<’•'*> =i-J 
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where R: is the radius of the Moon and r is 
the variable radius 


erfc X = 1 — erfx; erfx = 


y/irj 


The bracket in (1) is the first term of the 
series of the determinate Green’s function 
(1.6) (see Appendix I) : 


r (2n + 1) R - r ~ 


- er„ r(2» + .l) 

L V4fc (( - t) J 


The first term in the expression for tempera- 
ture gives the temperature of the central re- 
gion and the second gives the region of heat 
outflow to the surface. The region of outflow 
begins at depth: 

D = R-r<3\/ikt (2) 

If there is no convection, the region of out- 
flow, D, for Earth, Mars, and Venus does not 
reach the center, while D for the Moon 
reaches the center. 


Analysis of uniform source models of the 
Moon, with a variable initial temperature 
To (r) , is useful in searching for the required 
models of lunar accretion that would provide 
a high early temperature of the upper layers 
of the Moon. This analysis should give the 
required age of the mare basalts and should 
provide rapid melting and rapid differentia- 
tion of the upper few hundred kilometers and 
a relatively cold lithosphere at the present 
time. The mineralogy of the lunar rocks in- 
dicates a phase of high temperature and low 
pressure as suggested by Ringwood (ref. 26) 
in a semiqualitative conclusion (even before 
the Apollo missions) that the initial tempera- 
ture curve had maximum values in the upper 
200 to 300 km and a gradual drop toward 
the center (figs. 1, 2). This curve was used 
in the calculation of Kaula (ref. 29) during 
interpretation of the lunar gravitational 
field. Ringwood (ref. 30) intuitively pre- 
dicted a curve that is close to the curve for 
To(r) given by Hanks and Anderson (ref. 
31). 

Homogeneous models of the Moon, with 
such nonuniform initial curves with maxima. 


T^'K 



Fi^re 1 . — Variants of initial temperature distribu- 
tions in the primitive moon: 1 — from Ringwood 
(ref. SO), 2 — from Mizutani (ref. 27), S — from 
Hanks and Anderson (ref. SI), i — MacDonald 
(ref. 9), Ruskol (ref. S) (estimated as the result 
of collision of protomoons), 5 — Urey and Mac- 
Donald (ref. 28). 



Figure 2 . — Effect of Resorption of Initial Heat with 
Passage of Time in the Lunar Interior in Absence 
of Convection (curve 1, fig. 1 is taken as To). 
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indicate a possible thermal evolution that is 
dependent on the influence of initial pro- 
cesses. In particular, they permit understand- 
ing of the slowness of the disappearance of 
the maximum To{r) and the slowness of 
propagation of the initial heat inward to the 
planet. The latter conclusion follows from ex- 
amination of equation (2), which gives the 
rate of penetration of temperature by means 
of heat conductivity in the absence of convec- 
tion. Initially, the interior of the Moon 
apparently was cold. The task consists of as- 
certaining the manner in which the tempera- 
ture of the central regions of the Moon below 
1000 km could reach at the present moment 
high values sufficient for formation of the 
aseismic “asthenosphere,” the level at which 
deep moonquakes stop, according to data 
from low-frequency seismometers (ref. 18). 


Evolution of A Heterogeneous 
Moon 

Having selected a uranium content corre- 
sponding to the observed heat flux, U = 0.05 
ppm in the lithosphere and 0.5 ppm in the 
crust, we obtain the following models of the 
lunar thermal history (fig. 3). The tempera- 
ture curves are compared with the solidus 
curve of mare basalts, as determined by Ring- 
wood and Essene (ref. 26), in the zone from 
100 to 700 km, assuming no convection. The 
lithosphere thickens for 1.2 b.y. at a rate of 
160 km/b.y. Melting proceeds from shallow 
depths to great depths during the filling of 
the mare (3.4-3.1)-(3.8-3.6 AE), in agree- 
ment with the depth of origin for mare ba- 
salts, according to petrogenetical theory (ref. 
32) . Moreover, a sufficiently thick lithosphere 
is provided as required to support the stresses 
connected with the mascon gravitational 
anomalies. Between 1 and 2 AE there oc- 
curred conditions for the existance of a 
convecting core of undifferentiated lunar 
material which contained approximately the 
original concentrations of radioactive ele- 
ments. The core grew to a radius of 1200 km 
at 3 AE, then slowly decreased because the 


concentration of heat sources decreased with 
time. The chronology, mineralogy, and petro- 
graphy of the lunar mare basalts indicate a 
definite sequence of events in differentiation 
and fractionation of the material. For ex- 
ample, the mare basalts evidently had to be 
formed during the second thermal cycle, if 
there was one, i.e., 3.4 to 3.7 b.y. ago, at 
depths on the order of 200 to 300 km as a 
result only of endogenous factors. Crater for- 
mation apparently mainly ended during the 
first, most intensive, phase of differentiation 
and fractional crystallization. 

The heterogeneity of the internal structure 
of the current Moon is demonstrated by the 
low-frequency seismic experiment, data on 
the electrical conductivity distribution in the 
Moon, the existence of an ancient magnetic 
field, and gravitational data. A schematic 
section of the heterogeneous Moon was given 
by Anderson (ref. 33). The Moon consists of 
a crust, lithosphere, asthenosphere, and, pos- 
sibly, a core. The upper layer of the crust is 
composed of mare basalts. The lower layer of 
the crust is composed of gabbroic anortho- 
sites; under them is a layer of refractory 
crystals, and even plagioclase or pyroxenes, 
at depths of 60 to 80 km (ref. 34). The re- 
quirement for the absence of large-scale 
melting in the upper 1000 km imposes sig- 
nificant limitations on the thermal evolution 
of the Moon (ref. 6). 

The chronology of volcanic activity and 
the concentration of radioactive elements to- 
ward the surface are evidence of intensive 
differentiation at an early stage of lunar de- 
velopment (ref. 35). There are two points of 
view on the differentiation process. Differen- 
tiation may be the result of a relatively rapid 
accretion and subsequent general melting or 
of initial chemical stratification during the ac- 
cretion process. 

If the measured values of the heat flux 
0.74-0.67 X 10 ® cal/cm-sec (ref. 36) rep- 
resent the average heat loss characteristic 
of the entire Moon, the average uranium con- 
centration in the lunar lithosphere must be 
in the 50- to 80-ppb range. 

Attenuation of seismic waves depends 
strongly on temperature, showing a rapid 
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Figrure 3 . — Thermal evolution model for a heterogeneous Moon (this work) K/U-2000, Th/U 4, U — 
0.5 ppm in crust and U = 0.05 ppm in lithosphere. Initial temperature taken from Ringwood (ref. SO). 


increase with temperature and a sharp in- 
crease with the onset of melting (refs. 23 
and 37). Partial melting, with several per- 
cent of melt, can explain existing seismic 
data (ref. 23). The foci of moonquakes for 
which the depths were determined (18 cases) 
are concentrated in the 800- to 1000-km 
range; i.e., there is a concentration of lunar 
seismic activity in a zone 200 km thick. 

The silicate interior of the Moon has a 
solidus temperature from 1620 K to 1950 K at 
depths from 100 km to the center. This is 
compared with some thermal models (ref. 23) . 

Geological and geochemical processes on 
the Moon are considerably more primitive 
than the processes of crustal formation on 
the Earth. Magmatic activity on the Moon 
leads to eruption of magma and formation of 
gabbro-basalt-type rocks similar to the mag- 
matic rocks of the Earth. Differences are con- 
nected with the enrichment of the Moon with 
lithophile, refractory element, and the loss 


of a part of the siderophile elements. The 
depths of formation of magma chambers on 
the Moon, the rate of rise of magma, the 
temperature, and the diversity of surface 
crystallization conditions are reflected in the 
structure and texture of the rock (ref. 38). 
High concentrations of FeO, Ti02, and other 
elements reduce the crystallization tempera- 
ture of the magma. The ratios of close iso- 
morphic pairs of chemical elements K/Rb, 
Ca/Sr, Th/U approach the same ratios as in 
the tholeiite basalts of the Earth (refs. 22 
and 39). 

One of the hypotheses of the origin and 
evolution of the heterogeneity of the Moon 
which deserves attention is the hypothesis of 
convection (ref. 40). The following expres- 
sion for liquid flow of a cellular structure is 
obtained between the Nusselt (Nu) and Ray- 
leigh (Ra) criteria which describe convec- 
tive heat transfer: 
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Nu = 0.205 

Nu =4-1 Ra = -^ Al*; (3) 

a is the heat transfer coefficient, I is the 
thickness of the liquid layer, \ is the thermal 
conductivity coefficient, a is the thermal 
diffusivity coefficient, v is the kinematic vis- 
cosity coefficient, p is the bulk expansion co- 
efficient, g is the acceleration of gravity, and 
A is the temperature gradient in the layer. 

Experiments with various liquids at 4000 
< Ra < 10® have shown that equation <3) 
describes the test data well. The value of 
Racrit (which determines the start of the con- 
vective process) proved to be 3100 according 
to the experimental data. Thus, the convec- 
tive process starts at the moment when 
I > ^crit, where 

l„a = 4 / 3100 ^ ( 4 ) 

\ g p A 

Let the wall temperature be Ti and T 2 and 
(Ti < T 2 ). The magnitude of the heat flux 
from the hotter wall will be 

Q = a (Ti - T 2 ) (5) 

On the other hand, 

Q = A./ X A (6) 

The quantity A can be assumed to be equal to 
the melting temperature gradient, since it is 
known from geochemistry that moving melts 
are not superheated. Besides, this follows 
also from physical considerations, since the 
values of Ti and T 2 are essentially the melt- 
ing temperatures of rocks. 

Thus, we write relationship (3) in the 
form 

Q = lAa (7) 

Then, from relationships (3, 6, and 7) we 
obtain 

A.//A = 0.205 RaV-* (8) 

The values of Ae//A were calculated for a melt 
under Earth, Mars, and Moon conditions. 

Parameters adopted in the CGS system are 
shown in Table 1. 

Parameters of Thermal Models 

The basic parameters of thermal models 
are the coefficients of heat transfer, heat ca- 


Table 1 . — Parameters for Calculation by 
Equation (8) 


g, 10* 

a - 0.01; A - 3x10“ 

p, 10“ 

^crlt f 
V - 10“ 

(km) 

V - 10“ 

Earth 

1 

2.0 

47 

84 

Moon 

0.16 

4.0 

64 

113 

Mars 

0.39 

3.5 

52 

93 


pacity, latent heat of phase transitions, and 
the heat sources. The coefficient of heat trans- 
fer includes the coefficient of molecular 
(phonon) heat conductivity of a solid A</>, ra- 
diant heat transfer Ar, as well as the effective 
heat transfer by convection A*/. A family of 
and \4> + Ar curves for the presumed rocks 
of the Moon are presented in figure 4. At high 
temperatures, the coefficients change between 
the bounds of 0.003 and 0.012 cal/cm • s • K, 
and it is very likely that, deep in the interior 
of the lunar lithosphere, the value of the 
thermal conductivity is close to constant, at 
0.07 cal/cm • s • K = 2.8 W • m-‘ • K'^. 

The unique thermal properties of the low- 
conduction surface layer of the Moon must 
play a significant role in the heat balance 


T. C 



Figure 4 . — Thermal conductivity coefficients for lu- 
nar rock (curves 1-6 correspond to Toksoz and 
Solomon (ref. 2S); curve L, this work). 
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Figure 5. — Model of step toise change in thermal conductivity used in calculations of evolution of heterogen- 
eous Moon. 


and evolution of the heat flux. In connection 
with this, we have studied in greater detail 
than has been done previously the thermal 
conditions of contemporary models of the 
Moon, with the sharp change in thermal con- 
ductivity taken into account by means of the 
maximum possible decrease in the interval of 
temperature calculations close to the surface. 
Our initial model of thermal conductivity 
distribution is given in flgure 5 and table 2; 

The thermal isolation layer at the surface 


Table 2.-Initial Model of Thermal Conduc- 
tivity Distribution 


Depth, km 

Thermal Conductivity X, 
cal/cm • s • deg. 

0-0.250 

7 X 10-‘ 

0.250-0.760 

7 X 10-‘ 

0.750-1.000 

7 X 10-’ 

Ikm-lOkm 

7 X 10-“ 

lOkm-lOOOkm 

7 X 10-' 

1000km-1740km 

7 X 10“ 


of the Moon is of significant value for pro- 
longed retention of the melt in the interior 
of the Moon. If A = const = 0.007 cal/cm • 
s • deg., it persists for 300 m.y. ; if X is accord- 
ing to table 2, it persists for 1.3 AE. 

The model of heat generation in the in- 
terior of the Moon is determined by the con- 
tent of uranium, thorium, and potassium. The 
Luna and Apollo missions demonstrated that 
the chondrite model of the Moon is not cor- 
rect (ref. 22). Several variants of U, Th, and 
K content were tried in numerical calcula- 
tions. 

For the case of uniform distribution of 
heat sources, variants with the maximum B 
and minimum A in content were taken: 

A: U = 0.5 ppm 
B: U = 0.8 ppm 

After formation of a crust with layers 20 and 
60 km thick and a lithosphere (1000 km), 
redistribution of sources 'leads to approxi- 
mately a twentyfold depletion of uranium, 
according to the model of Vinogradov (ref. 
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22). According to these data, the concentra- 
tions of the elements in the upper 20 km was : 
U = 2.3 X 10-" g/g ; Th = 0.9 X 10“^ g/g ; 
K = 0.83 X 10-*%. 

After stratification, the following remained 
in the lithosphere: 

U = 0.43 X 10-* g/g; Th = 0.38 X 10~^ 
g/g; K = 0.83 X 10-*%. 

During melting, the physical parameters of 
the phase transition boundaries are dis- 
rupted. This disturbs the stability of the cal- 
culation schemes, if special smoothing 
procedures are not used. The effects of melt- 
ing and convection on temperature were 
modeled using a technique described by Bu- 
dak et al. A smoothing function of the physi- 
cal parameters at the phase transition 
boundary was introduced. When the tempera- 
ture T“ (r,t) becomes higher than the melt- 
ing temperature, the thermal conductivity 
coefficient is replaced by an effective heat 
transfer coefficient kef = 0.1-0.05 cal/cm • s • 
° C, symbolizing convection. 

The heat sources are the most important in 
calculations of the thermal history. The 
abundance of isotopes measured in lunar 
samples and by means of orbital gamma-ray 
observations (ref. 41) leads to the conclusion 
that the concentrations of uranium and tho- 
rium in the surface layers of the Moon are 
increased over those in the deep layers and 
are generally higher than on Earth. At the 
same time, the Moon is depleted in volatiles, 
including potassium. While for Earth K/U = 
10 000 and for chondrites K/U = 80 000, we 
have only K/U = 2000 for the Moon. Using 
the ratio K/U = 2000, we obtain the aver- 
age generation of heat H (r,t) in the interior 
of the Moon as a simple function of the bulk 
uranium content U (to) : 

H (r,t) = p U (to) 1.35 [0.73 + 0.20 Th/ 
U + + 0.25 X 10-" K/U] 
p = 3.34 g/cm® 

Uo = 7.3 • 10-« g/g = 7.3 ppm 
After stratification during the first billion 
years, the crust of the Moon consists of 20 
km of basalts and 40 km of andesites; the 
Moon is considered to be three-layered, ac- 
cording to Anderson (ref. 33). 

The average content of uranium is between 


those of howardites and eucrites and deter- 
mines the average value of the global heat 
flux which is 33 erg/cm* • s (ref. 42) . The de- 
cay constants, decay energy, and isotope 
abundances are taken in accordance with ref- 
erence 43. We disregard the value of tidal 
friction, including it in To (r ) . 


Direct and Indirect Determination 
of Heat Flux and Structure of 
the Lunar Crust 


Before the flights to the Moon of the Luna 
and Apollo 15, 16, and 17 spacecraft, the heat 
flux on the Moon was determined by study of 
its natural radio emission and extraction of 
the constant components that determine the 
temperature gradient in the interior. Pre- 
cision measurements of the constant compo- 
nent of the radio temperature, at two 
wavelengths, pi and p 2 , perpendicular to the 
surface, give an expression for the tempera- 
ture gradient (refs. 15 and 44). 


grad T(z) = 


dTeo 1 

dp (1 — ^) mlt 


where T is the average effective temperature 
over the disk, depending on p; is the al- 
bedo; and le/lt is the ratio of the depth of 
penetration of electromagnetic and thermal 
waves, Si/p = m. 

The quantity U, the depth of penetration 
of electromagnetic waves, is proportional to 
wavelength p and inversely proportional to 
tan A, the tangent of the loss angle (ref. 44) : 


le = p/(2,ir Ve tanA) = d p 
where t is the dielectric constant and d is a 
proportionality coefficient, depending on the 
properties of the substance in which 
tan A = 2<r/tf 

Here, a is the effective electrical conductivity 
at a given frequency /. Multiplying grad T 
by the thermal conductivity coefficient A, we 
obtain the heat flux. 


_ dTeo Vtt 

dp (1 — Ej) my \/t 
Using the model of heterogeneous lunar 
surface structure (rocks with pi covered by 
a porous layer with po), Tikhonova and 
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Troitskiy (ref. 15) predicted the relatively 
high value of the heat flux from the interior 
of the Moon: 0.7 X 10~« < <7(j < 0.95 X 10-® 
cal/cm^ • s. 

Direct measurements of the heat flux in 
boreholes provide significant information on 
the heat balance and evolution of the Moon. 
Direct measurements of the heat flux from 
the interior of the Moon, with and without 
introduction of regional corrections (Apollos 
15 and 17) , showed that the value of the heat 
flux is at the upper limit of the expected 
values given on the basis of geochemical 
models, the intensity of radio emission, and 
calculations of the thermal history of the 
Moon. The two groups of measurements in- 
dicated were made at the margins of Mare 
Serenitatis and Mare Imbrium, at the boun- 
dary of the assumed mascon basins (fig. 6). 
The thick lunar crust contains significant 
irregularities. Particularly large irregulari- 


ties are concentrated on the visible side of the 
Moon. 

Similar boundaries are connected with 
change in the physical-chemical properties. 
In this connection, it is natural to expect a 
change in the thermophysical properties in 
the transition from the mare regions to the 
highlands. In fact, investigations by Luna 21 
(Lunkhod 2) in the region of the eastern mar- 
gin of Mare Serenitatis (as well as Lunokhod 
1 in Mare Imbrium) showed that, according 
to the data of the RIFMA-M instrument, the 
iron content in the highlands is less than in 
the mare. According to radar observations, 
the region of measurement of the thermal 
flux was close to a “black spot region.” The 
blue tint of these spots also is evidence of a 
high concentration of iron and titanium ions, 
with which variations in the coefficient of 
thermal conductivity and refraction of the 
heat ffux may be connected. Terrestrial ex- 
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perience with heat flux measurements indi- 
cates the unreliability of judgments as to the 
average heat loss, q, on the basis of one or 
two measurements. 

On the other hand, the absence of even a 
small region of water and appreciable tec- 
tonic activity permits one to assume that 
local variations in the heat flux on the Moon 
will be less than on Earth. However, the 
specific factors of disintegration and erosion 
of the surface of lunar rock (meteorite im- 
pacts, thermoelastic stresses as a result of 
fluctuations of the surface temperature, the 
effect of the mascons as inclusions of extrane- 
ous thermal conductivity) could be reflected 
in the history of the heat flux. The phenom- 
ena of disintegration of the lunar surface 
rocks by meteorite impacts should have intro- 
duced local changes in the energy balance 
of the surface. 

The places of measurement of the heat flux 
in the boreholes on the Moon are presented 
in figure 6. They are located at the edges of 
Mare Serenitatis and Mare Imbrium. The 
Apollo 15 borehole is in the region of the 
Appennine Mountains, between Mare Sereni- 
tatis and Mare Imbrium and distant from 
three craters, Autolycus, Aristillus, and Ar- 
chimedes. The landing site of Apollo 17 is 
located in the Taurus-Littrow Valley. There 



Figrure 8 . — Change in relative value of heat flux 
q/q„ above horizontal conducting inclusion (A) 
(cross hatched) with thermal conductivity three 
times the thermal conductivity Xi of the surround- 
ing rock, and above a thermal insulating inclusion 
(B). Curves 1, 2, and S correspond to depths of 
burial: (1) H = 0, (2) H = i km„ (S) H = 2 km. 


is a cross section of the latter, according to 
seismic sounding data, given in figure 7 (ref. 
45). We made an attempt to determine the 
shape of the heat flux anomaly at the edges 
of the mascon, on the basis of numerical mod- 
eling of a two-dimensional structure, con- 
taining inclusions with high or low thermal 
conductivity as functions of varying degrees 
of burial (fig. 8). Data on thermal conduc- 
tivity of the lunar rock and observed varia- 
tions in them were used. Two mascon models 
were used (figs. 8 and 9), with and without 
allowance for the relief. In accordance with 
conclusions from the data of Luna 10, it was 
assumed that the surrounding surface rock 
of the Moon was basaltic. 

The heat-flux anomaly has the appearance 
of sign-changing curves on the boundaries 
of trough-shaped (fig. 9) and rectangular 
(fig. 8) irregularities. With movement of an 
inclusion to the surface, a false anomaly at 
its edge can exceed 20 percent, with a three- 
fold thermal conductivity drop. The size of 
the anomaly can increase due to a sudden in- 
crease in uranium concentration. On the 
lunar highlands and mare, Langseth et al. 
(refs. 16 and 17) used the observed q and 
more accurately determined the structure of 
the cover. 

The absence of water and of noticeable 


Fi^re 7 . — Cross section of the structure of the lunar 
crust in Apollo 17 landing region, from seismic 
data; line A-B, assumed fracture (ref. U5). 
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Figure 9 . — The thermal effect of a conducting mascon 
(Xi = iO \ij, buried 20 km in a poorly conducting 
layer (\i = Xi/250) and with a high temperature 
Ti at the level of the compensating layer. 

tectonic activity on the surface of the Moon 
allows one to assume that local variations in 
the heat flux on the Moon will be less than 
on Earth. However, disintegration and ero- 
sion of the surface layers by meteoritic im- 
pacts has severely disrupted the surface 
layer. Whether or not the direct measure- 
ments of the heat flux (refs. 16 and 17) by 
Apollo 15, 16, and 17 are representative in- 
volves data on the structure of the upper 
layer. 

Research of Troitskiy et al., (ref. 44), 
Tikhonova and Troitskiy (ref. 15), and 
Langseth and Keihm (refs. 16 and 17) was 
devoted to study of the nature of the mate- 
rial and its density and heat conductivity. 
Signiflcant change in the heat conductivity 
of the regolith and deep underlying layers 
was demonstrated in these researches by in- 
dependent methods. In the work of Langseth 
and Keihm (ref. 36), the condition of a 15- 
cm layer of the regolith is studied in detail. 
At the same time, in work on the thermal 


history of the Moon, one studies the distri- 
bution of the thermal parameters on a 
considerably rougher scale which is unsatis- 
factory for calculation of the heat flux. 

Calculation Procedure 

The analytical solution for temperature 
and heat flux presented in Appendix I is a 
general solution that is universal for spheri- 
cal planets of any radius, for stages when 
the distribution of heat sources and thermo- 
physical parameters can be assumed to be 
uniform but, permitting change with time 
(e.g., a decrease in the concentration of ra- 
dioactive elements because of their contin- 
uous decay). 

In all later stages of development of 
planetary bodies, with the exception perhaps 
of the asteroids, the distribution of properties 
and of heat sources in particular is irregular, 
and differentiation, melting, and convection 
are included. The problem is solved with 
computers, using the method of finite differ- 
ences. Any temperature distribution given by 
a computer is a partial solution. In principle, 
there can be an infinite number of such so- 
lutions from the point of view of thermal 
conductivity theory. In our numerical exam- 
ples, the thermal conductivity equation was 
approximated by a finite difference scheme, 
with steps h along the axis r and along the 
time axis t (Appendix II). During melting, 
part of the thermal energy must have been 
absorbed as the latent heat of fusion. 

Thermal conductivity problems, with phase 
transitions on the moving interface of two 
phases taken into consideration, are known 
in thermophysics as the “Stefan problem.” 
Current models of the thermal evolution of 
the surface molten layer of the Moon (refs. 
23 and 13) show that this layer is gradually 
displaced inward, with melting occurring on 
the lower boundary, while solidification and 
crystallization occur on the upper. In this 
manner, we must consider the Stefan prob- 
lem for two moving boundaries: the upper 
z = Xi (t) and the lower z = X 2 (t). It is 
controlled by the heat balance conditions in 
the form 
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dXi(t) _ dTi 


dt 


— Ai ■ 


dXj{t) _ 


dt 


— — X2 


dz 

dTi 

dz 


— X 2 


+ X 3 


dT, 


dz 

dTs 

dz 


z = x 2 


«= *2 


where z is the depth ; Ai and X 3 are the ther- 
mal conductivities above and below the mol- 
ten layer ; X 2 is the coefficient of heat transfer 
or effective thermal conductivity inside the 
molten layer, possibly with inclusion of con- 
vective transport; L is the latent heat of fu- 
sion ; and Ti is the temperature in the upper 
solid layer, satisfying the boundary condi- 
tions: 

TAz)\ =0 

U = 0 

Tr{z)\ =T„\_X,{m 

\ Z - Xl (t) 

Tm = To + y 2 Z 

Here, Tm is the melting temperature, and 72 
is the melting temperature gradient. Let T 2 
be the temperature inside the molten layer, 
satisfying the conditions: 

T2(z,t)\^^^^^ T„ [xi (f)] 
T2(Z,t)\^^^^= [X2 (f)] 

T 3 is the temperature in the lower layer be- 
low the level of the melt, satisfying the con- 
ditions : 

T3(z,t)\ = [X 2 (f)] 

dTs 

dz , = ,2 “ ■>'2 

The initial conditions for the moving boun- 
dary are the following : 

The Stefan problem is nonlinear. The non- 
linearity increases if the dependence of ther- 
mal conductivity on temperature is taken 
into consideration. In order to carry out auto- 
matic soothing of jumps in the physical pa- 
rameters on the phase transition boundaries, 
we use the smoothing method of Budak et al. 
and we introduce a change of the variables 

T ^ U + T* 

and the specific heat content function (see 
Appendix III) 


H (u,r) = S lcU)pU) +.LpSU - r*)]d^ 


The difference analog of the thermal con- 
ductivity equation is given in Appendix III. 
The value of the latent heat is assumed to be 
L = 400 J/g. The surface temperature is as- 
sumed to be constant during the life of the 
Moon and equal to —20° C. 

The temperature distribution should be 
satisfied by the data of a comparatively high 
heat flux (refs. 16 and 17). The value of q, 
previously calculated from theory for the 
present and with the age of the Moon being 
4.6 X 10® yr, is 0.36 X 10~® cal/cm^ • s (ref. 
2) . A molten layer forms at 2.1 X 10® yr, at 
a depth of 500 km, with a 150-km width and, 
at 3.2 X 10® yr, the melting level reaches 300 
km and is maintained up to the present time. 

The hot state of the Moon could be ex- 
plained in this manner, but the model does 
not completely satisfy the heat flux and the 
heterogeneous seismic structure. Better 
agreement with the seismic model and geo- 
chemical evolution is obtained by assuming 
higher uranium content up to a value of 60 
ppb in a homogeneous Moon, which corre- 
sponds to the models of Toksoz et al. (ref. 46) 
and Toksoz and Solomon (ref. 23). In figure 
3, with an initial formation of a melt, the 
effect of the latent heat on the phase tran- 
sition boundaries was taken into account: 
absorption of 400 cal/g at the upper bound- 
ary and emission on the lower boundary of 
the melt, in accordance with Stefan’s law for 
a moving boundaiy. The latent heat was in- 
corporated in the thermal conductivity equa- 
tion through the specific heat content by the 
method of Budak, and the solution of Stefan’s 
problem then has the form of Appendix III. 
At the onset of melting, the heat transfer in- 
creases, and the coefficient A increases ten- 
fold. If one assumes differentiation at the 
very Start and a crustal thickness of 60 km 
with uranium concentration 0'4 X 10~® g/g, 
then a melt forms in the lithosphere at a 
depth of 100 - to 600-km and quickly expands, 
and the greatest heat flux reaches 1.17 X 
10“® cal/cm^ • s (1 AE) and is 0.6 X 10~® 
cal/cm® • s at present, which is in agreement 
with the heat flux predicted by Tikhonova 
and Troitskiy (ref. 15) and measured by 
Langseth and Keihm (refs, 16, 17, and 36). 
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Consequences of the Thermal 
History of the Moon 

study of moonquakes and lunar tectonics 
by high-frequency seisnaometers has shown 
that thousands of very small lunar signals 
recorded by the Apollo 14, 15, and 16 stations 
are correlated with meteorite impacts, as 
well as with the sunrise and sunset periods. 
This correlation of micromoonquakes with 
the thermal cycle on the surface indicates 
that a considerable fraction of the micro- 
moonquakes are of a thermal nature (ref. 
47) . There are two likely mechanisms of such 
periodic thermal moonquakes: (a) fractur- 
ing of the lunar rock and (b) slippage of 
lunar soil along slopes, caused by the periodic 
thermal stresses. The moonquake signals re- 
corded by high-frequency seismometers are 
completely different from locally generated 
moonquakes. In conformance with current 
seismic models of the structure of the Moon, 
deep moonquakes arise at the base of a thick 
elastic shell (lithosphere) , immediately above 
the boundary of the relatively weak central 
zone of the asthenosphere. Surface tectonic 
characteristics indicate the possibility of a 
slight tension on the lithosphere of the Moon 
in the past when the interior of the Moon 
became warmer. However, the absence of 
moonquakes at shallow and intermediate 
depths at present leads to the assumption 
that no appreciable compression or dilation 
rates exist now. Therefore, the Moon should 
be close to thermal equilibrium, and the rate 
of heat loss by radiation of heat from its 
surface should be approximately compensated 
by the rate of internal heat generation. 

Seismic and tectonic activity of the Moon 
is extremely low compared with Earth. The 
Moon is characterized by an inactive outer 
layer 1000 km thick. The seismic hypothesis 
can be tested from the point of view of ther- 
mal processes by calculation of the thermo- 
elastic state of the Moon, with allowance for 
the fact that its upper layer has already 
passed through a molten state and that the 
central layers are now in a softened state 
close to partial melting or the beginning of 
melting, i.e., the solidus. 


The great thickness of the lunar litho- 
sphere, compared with that of Earth, most 
likely explains the strongly differing tectonics 
of these planets. 

An estimate of the propagation of elastic 
deformations and stresses that result from 
heating and cooling of the Moon can be car- 
ried out for models of an elastic, spherical 
layer placed on the surface of a molten layer. 

For a maximum estimate of the expansion 
of the perimeter of the Moon, the elastic 
model can be used. The change in radius of 
an elastic Moon during its thermal history 
can be estimated by the formula 

dEij a f . HRaA 

where q is the heat flux at the surface, a is 
the coefficient of thermal expansion, and H 
is the heat generation. In 1960, the rate of 
expansion of the Moon was estimated by 
MacDonald, and in 1968 by Lyubimova, on 
the basis of chondrite models of the Moon 
and of a comparatively low heat flux taken 
from calculations of the thermal history. We 
have available at present independent, di- 
rectly measured values of and U. We 
substitute the following values of the flux 
% = 28 erg/cm^ • s, according to Langseth 
et al. (refs. 16 and 17) and heat generation 
based on the average uranium content U = 
60 ppm and the value a = 3 X 10”® deg”^ 
(ref. 23). By substituting these numerical 
values, we obtain 10 ® cm/yr for the maxi- 
mum rate of change in radius of the Moon 
in the last billion years of its evolution. The 
rate of change in the radius of the Moon 
was greater in the past than now. Deforma- 
tion u and tangential stresses and t on the 
surface r — R, in the presence of a molten 
layer, can be determined from the equation 
of elastic equilibrium in a spherical layer 
bounded on the inside (r = a) and outside 
(r = R) by spherical boundaries. The fol- 
lowing are conditions on these boundaries: 
the upper one is free, and a normal stress 
is assumed on the lower one, designated con- 
ventionally by y. For a given temperature 
distribution, we obtain deformation 


342 


COSMOCHEMISTRY OF THE MOON AND PLANETS 


and also, 

dr r 


--^(1 + ,) (1-25) 
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i.e., the increase in radius m (72) and tangen- 
tial stress t(72) on the surface 
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The formation of fractures on the surface 
of the Moon is possible only upon reaching 
T, the limit of the strength of the material. 
From the relationship between u and f just 
obtained, the magnitude of the fractures can 
be estimated, if the value of the critical stess 
T„it is substituted in the left side. We assume 
that on the boundary r = a <!„ = y. Since 
the radial component of the stress in a liquid 
is given by the formula 

Vo = Vt = 0 

on the boundary, r = a and = y 
7 = p (a) =2 - -g- divlf ^ 

where pi is the density of the melt. The for- 
mula for determination of the increase in 
perimeter of the Moon will have the form: 
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Calculations show that the perimeter of the 
Moon increased during the first half billion 
years and then began to decrease until the 
time of formation of the asthenosphere. If 
the assumption of the existence of the strong 
convective heat transfer in the region of 
the partially molten asthenosphere of the 
Moon is introduced, an increase in the perim- 
eter again resulted 2 to 3 b.y. ago. It is as 
though the molten layer expands the Moon 
from inside, and the increase in the perim- 
eter ensures formation of crustal fractures 
through which the magma erupts. 


Conclusions 

The principal attention has been given in 
this work to an effort to relate the values of 
the observed heat fluxes to the thermal evo- 
lution of the Moon, considering the fine 
structure of the outer crust which has the 
properties of a uniquely insulating material. 

In calculations of the history of the heat 
flux, extremely small depths (steps in the 
calculation network of 250 m in all) were 
introduced into analysis of the numerical 
determinations of temperature differences 
(on the order of hundreds of meters), with 
allowance for abrupt differentiation in the 
thermophysical properties and for heat gen- 
eration. 

In this work, initial melting of the upper 
250 km was assumed, with allowance for 
abrupt stepwise differentiation of the ther- 
mophysical properties and heat generation, 
as well as for the rapid fractionation of mat- 
ter. The analysis included the effects of 
absorption and emission of heat at the bound- 
aries of a moving, molten layer, with effec- 
tive internal convection and rapid transport 
of the excess heat from the upper boundary 
to the lower one. 
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A new technique based on the Stefan equa- 
tion is introduced for calculation of the ab- 
sorption and emission of heat that results 
from the effect of the latent heat of fusion 
of phase transitions occurring at two moving 
boundaries of a melt. 

From the analysis presented in this paper, 
the following preliminary conclusions are 
made: 

1. The introduction of thermal conduc- 
tivity, decreasing toward the surface, 
increases the thermal insulation effect 
of the upper cover of the Moon and 
facilitates retention of the melt in the 
interior of the Moon for a period of 
1.3 b.y., rather than 300 m.y. as for a 
uniform thermal conductivity of 0.007 
cal/cm • s • °C. 

2. A uniform distribution of uranium with 
a concentration of 60 ppm does not 
reflect the observed values of the heat 
flux because y-spectrometry data indi- 
cate a sharp concentration of radioac- 
tive elements in the crust and their 
enrichment in the lunar highlands. 

3. A distorting effect of the mare basin 
mascons on the observed heat fluxes is 
not excluded due to the different nature 
of these basins whose substance is more 
heat conducting and contains increased 
percentages of iron and titanium. The 
greatest distortions are expected near 
the edge of such basins where the heat- 
flux measurements were made. It is 
presently difficult to accurately predict 
the magnitude of the expected anom- 
alies because of the shortage of data 
on the thermal conductivity coefficient 
of the matter making up the mascons. 

4. The thermoelastic history of the Moon, 
analyzed with allowance for the ex- 
panding effect of the molten layer be- 
low the lithosphere, shows that the 
increase in the Moon’s perimeter was 
sufficient for the formation of deep 
fractures through which lavas erupted 
and filled the mare basins. The erup- 
tion of lava was stimulated in the early 
history of the Moon by bombardment 
of the surface with meteorites. 


5. During the life of the Moon, the emis- 
sion of heat into space has gradually 
declined. At 1 AE after formation of 
the Moon, the heat flux reached 1.17 X 
10-® cal/cm* • s ; after 3 AE, it decreased 
to 0.64 X 10"® cal/cm* • s and almost 
stabilized. Since then subsequent de- 
crease has not exceeded 10 to 15 per- 
cent. 


Appendix 1: Analytical Solution 


Heat calculations are carried out using a 
nonstationary thermal conductivity equation 
for a spherically symmetric planetary body 
(Moon, Mars, Mercury, asteroids, etc.), con- 
taining internal heat sources 


Cpp 


dT ^ J ^ ■ 

dt r* dr 


r^(Kp + 


where Cp is the specific heat capacity, T is 
the temperature at point (r,t), p is the den- 
sity, and H(r,t) is a function of the heat 
source, which usually implies the long-lived, 
naturally radioactive isotopes U**®, U**®, 
Th***, and K. The energy of gravitational dif- 
ferentiation given off in the process of evolu- 
tion of the Earth, the latent heat of phase 
transitions, and other things can also be in- 
cluded in H(r,t) (ref. 2). 

The physically based boundary conditions 
are 


T(r = R,t) = Const \ 

^{r = Q,t)=Q ^ ( 1 . 2 ) 

T{r,t = 0) = To(r) ) 

For these conditions, the classical form of 
solution applicable to the radius of the Earth 
was obtained by Lowan (1936), Urey (ref. 
7) , and Kopal (ref. 48) . 

The expression for the temperature inside 
a sphere of radius R obtained by these re- 
searchers is presented in the form of an in- 
finite series by trigonometric functions and, 
as Van Orstrand showed (ref. 49), the series 
converges very slowly. The use of Green’s 
function, with application of the method of 
sequential reflections, gives considerably 
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more rapid convergence of the series. It also 
represents the temperature inside the sphere 
and the thermal flux on the surface in the 
form of a formula, permitting a clear physi- 
cal interpretation. Construction of a solution 
for temperature is carried out in the follow- 
ing manner. The substitution of variables is 
u = rXT{r,t) (1.3) 

The boundary conditions will then be 
u(z,0) = 0;u(R,t)= 0]u (0,t) =0 (1.4) 

The expression for u (r,t) is expressed 
through Green’s function G: 

t X 

u(r,t) ~ j j dr' dr (1.6) 

• • 

where G (r,t,r',T) is Green’s function for the 
segment [0,R'\ 

G{r,t,r'T)=% \_Go{r,t,2nR+r'X) — 

n= — 00 

Go(r,f,2-ni2— r'r)] (1.6) 

in which 

G.(n*r,r) = 2 vS7(^) 

(1.7) 

Further, the variables are replaced, 

r±2nR±r'=y (1-8) 

and the expression for the temperature in- 
side a spherical body takes the form 

I / os X-r 

* x-r 

00 

(r - 2nR + y)exp dy+ J 


R+r \ 

I ir + 2nR -y)exp[ - dy [ dr 

(2n+l> R-r / 
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After transformation of the sums and in- 
tegrals, we have 

t t 
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■erfc 
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The first term defines the temperature in 


dr 


the center; the second the magnitude of the 
outflow of the heat to the surface. Analysis 
of the convergence permits determination of 
the number of terms of the series n >m, 
starting with which, all subsequent terms of 
the series, with m > 3 \/Kt/R, can be dis- 
carded. 

With thermal diffusivity of the rock on the 
order of K = 0.01 cmVs, it holds true for 
the terrestrial planets and for intervals of 
time less than 5 X 10® yr that it is sufficient 
to take only one term of the series with n = 0, 
if the radius of the body is greater than 
R = 1700 km. 

The temperature in the center of the planet 
or asteroid of any radius and at any moment 
of time is given by the expression 

T(r = o,t) = lim T(r,t) = — / ff(r)dr - 

r-»0 Cp J 


where 


^ JL f-Mkl=r 

Cp \/tt J \/K{t—r) 

9 

r (2w+D® R^~\ 
L J 


£ 

dr 


H(t) = XHo exp( -\i(t-ro)) 


Appen(Jix 2: Numerical Solution 

Without Consideration 
of Melting 

The schemes commonly used previously 
had the shortcoming that, with small steps 
along the radius, too much machine time was 
required. Therefore, for example, MacDonald 
(ref. 8) had to limit the radius steps h to the 
order of 50 km; Toksoz and Solomon (ref. 
23) and Mayeva (ref. 50) to a value h = 20 
km ; and Ornatskaya (ref. 13) to about 5 km. 
In this case, the fine structure of the lunar 
crust and the unique nature of the low co- 
efficients of thermal conductivity of the outer 
crust of the Moon actually drop out of the 
analysis. Calculations of the surface heat flux 
at such h is extremely inaccurate. In this 
case, the effect of the irregular structure of 
the crust of the Moon is ignored. Limitations 
on stability are weakened by the use of im- 
plicit finite difference schemes. The error is 
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small for great depths, but it increases near 
the surface. In order to obtain a denser loca- 
tion of points near the surface, one can sub- 
stitute variables 

r = l-x2 (2.1) 

Then the thermal conductivity equation will 
be approximated by the scheme described for 
equal new variable steps. The calculation in- 
volves use of the standard trial-run method, 
a convenient method for solution of the re- 
sulting system of linear equations “linked” to 
each other. The finite difference analog of 
the thermal conductivity equation, with ir- 
regular distribution of heat sources H(r,t) 
and variable physical parameters p,A,Cp at 
each point of the calculation network, has 
the form 


° ciuMr,) I 


TJ 


Tj+Tj-i . 


Appendix 3: Consideration of the 
Latest Heat of Fusion. 


The difference analog of the thermal con- 
ductivity equation, containing a change in 
heat content due to absorption or emission 
of the latent heat of fusion at the phase 
transition boundaries, has the form in di- 
mensionsless variables 



+ LV.i)L 

where t' = t/L, r' - r/R 

hi = r'i - r'i-i 

, _ hi+hi^i 

hi 2 

^hi = l, K = L/R^ 
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